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PREFACE 


The Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing is a multiyear program of research, development, evaluation, and appli- 
cation of aerospace remote sensing for agricultural resources, which began in 
fiscal year 1980. This program is a cooperative effort of the U.S. Department . 
of Agriculture, the National Administration (U.S. Department of Commerce), the 
Agency for International Development (U.S. Department of State), and the 
U.S. Department of the Interior. 

The work which is the subject of this document was performed by the Earth 
Resources Applications Division, Space and Life Sciences Directorate, Lyndon B. 
Johnson Space Center, National Aeronautics and Space Administration and 
Lockheed Engineering' and Management Services Company, Inc. The tasks performed 
by Lockheed Engineering and Management Services Company, Inc., were 
accomplished under Contract NAS 9-15800. 
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1. INTRODUCTION 


The FRIEDMAN program is a FORTRAN-coded implementation of Friedman's nonpara- 
metric test for several related samples with one observation per treatment/ 
block combination, or as it is sometimes called, the two-way analysis of 
variance by ranks. This test is particularly useful when the assumptions of 
normal theory fail to hold for a randomized, complete block design. In fact, 
for this test, the only assumptions made are that the blocks are mutually 
independent and that within each block the data may be ranked from smallest to 
largest. A test statistic is then formed based upon the ranks, not the actual 
data. This statistic is then used to test the null hypothesis, as follows: 

Treatments have identical effects. 

H,: At least one treatment tends to yield larger 
observations than at least ono other treatment. 

For further oxtails, see reference 1. 

A detailed description of the FRIEDMAN program is given in this document. A 
test data set and its results are also presented to aid the potential users of 
this program. 
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2. PROGRAM DESCRIPTION 


2.1 GENERAL OVERVIEW 

The FRIEDMAN program calculates Friedman's T-statistic for data sets with a 
maximum of 500 blocks and 30 treatments, allowing for only one observation per 
treatment/block combination. This is done via the main calling program, MAIN, 
and five subroutines: INPUT, SORT, XRANK, TSTAT, and OUTPUT, The FRIEDMAN 

program utilizes an execution time format, providing for a great deal of 
flexibility in the configuration of the experimental data set. The program is 
written in FORTRAN IV-H and is currently implemented on the AS/3000 computer 
in the Earth Observations Data Laboratory, NASA Johnson Space Center. 


3. MAIN 


The MAIN program is the main routine of FRIEDMAN. This component reads from 
the terminal the number of card images per block, the number of treatments, 
and the execution time format. Subsequently, it calls the subroutines that 
actually calculate the T-statistic. 

CALLING PROCEDURE : Not applicable. 

INPUT PARAMETERS ; NC, NT, FMT 

OUTPUT PARAMETERS : M 

REFERENCED BY : Not applicable. 

SUBPROGRAMS REFERENCED ; INPUT, SORT, XRANK, TSTAT, OUTPUT 

INPUT/OUTPUT DEVICES ; Unit 16 (terminal). Unit 6 (disk, LRECL = 133, 

BLKSIZE = 133) 

LOCAL VARIABLES ; The following abbreviations are used in tables throughout 
this document. 

A = alphanumeric 
I = integer 
R = real 

Local variables for the MAIN program are as follows: 

Name Type Dimension Description 

XINPUT R 500 X 30 Input data set. 

RANK R 500 x 30 Array of ranked data. 

R R 1 X 30 Sum of ranks for each treatment. 
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Name 

Im. 

Dimension 

Descri ption 

DEV 

R 

1 X 30 

Array of squared deviations from 
the expected value of the treatment 
sum of ranks. 

WORK 

R 

30 

Working storage. 

NC 

I 

1 

Number of card images per block. 

NT 

1 

1 

Number of treatments. 

M 

I 

1 

Number of blocks. 

FMT 

A 

80 

Format of XINPUT. 
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3.1 INPUT 


The INPUT subroutine reads the exp\sri mental data set from the disk by reading 
in all treatments for a fixed block. Tliat is, treatments are arranged 
horizontally and blocks are arranged vertically. (See figure 1 for data 
configuration. ) 

CALLING PROCEDURE ; CALL INPUT(XINPUT, FMT,M,NT) 

INPUT PARAMETERS ; FMT, M, NT 
OUTPUT PARAMETERS ; X INPUT 

REFERENCED BY ! The INPUT subroutine is referenced by the driver routine MAIN. 
SUBPROGRAM REFERENCED ; None. 

INPUT/OUTPUT DEVICES ; Unit 5, the disk containing the experimental data set 
LOCAL VARIABLES; 


Name 

Jm. 

Dimension 

Description 

XINPUT 

R 

M X NT 

Input data set. 

M 

I 

1 

Number of blocks. 

NT 

I 

1 

Number of treatments. 

FMT 

A 

80 

Format of XINPUT. 
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Figure 1.- Input data configuration. 


3-4 


3.2 SORT 


The SORT subroutine sorts the treatment values for each block in ascending 
order. The algorithm used is a Shell sort (ref. 2). 

CALLING PROCEDURE ; CALL SORT(XINPUT,INDEX,I,M,NT) 

INPUT PARAMETERS ; XINPUT, I, M, NT 

OUTPUT PARAMETERS : XINPUT, INDEX 

REFERENCED BY ; The SORT subprogram is referenced by the driver routine MAIN. 
SUBPROGRAMS REFERENCED ; None. 

INPUT/OUTPUT DEVICES ; None. 

LOCAL VARIABLES ; 

Name Type Dimension Description 


XINPUT 

R 

M X NT 

On input, the raw experimental data; 
on output, the sorted experimental data. 

INDEX 

R 

NT 

Array recording treatment exchanges that 
occur during the sort. 

I 

I 

1 

Number of current block being sorted. 

M 

I 

1 

Number of blocks. 

NT 

I 

1 

Number of treatments. 

NDELTA 

I 

1 

Variable used to partition each block 
for sorting. 

NSWITC 

I 

1 

0 if no switch is made; 1 if a switch is mad( 

MAX 

I 

1 

NT - NDELTA. 

TEMP 

R 

1 

Temporary storage used when interchanging of 
treatment values occurs. 
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3.3 XRANK 


The XRANK subroutine, using the sorted experimental data, checks for ties and 
then assigns ranks. Using the array INDEX, XRANK positions treatment I in 
column I of array RANK. 

CALLING PROCEDURE ; CALL XRANK(XINPUT, RANK, INDEX, I, M, NT, WORK) 

INPUT PARAMETERS ; XINPUT, INDEX, I, M, NT 
OUTPUT PARAMETERS : RANK 

REFERENCED BY ;:The XRANK subroutine is referenced by the driver routine MAIN. 
SUBPROGRAM REFtRENCED : None. 

INPUI/O r/PUT DEVICES ; None. 


LOCAL VARIABLES ; 

Name Type Dimension Description 


XINPUT 

R 

M X NT 

Array of sorted experimental data. 

RANK 

R 

M X NT 

Array of ranks. 

INDEX 

R 

NT 

Array of treatment exchanges that occur during sort 

I 

I 

1 

Block number. 

M 

I 

1 

Number of blocks. 

NT 

I 

1 

Number of treatments. 

WORK 

R 

NT 

Working storage. 

NEQUAL 

I 

1 

Number of treatment values making up a tie. 

N FIRST 

I 

1 

Column number of first of a sequence of equal value 

LAST 

I 

1 

Column number of the last of a sequence of equal 
values. 
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3.4 TSTAT 


The TSTAT subroutine calculates Friedman's T-statistic based upon the contents 
of the array RANK. 

CALLING PROCEDURE : CALL TSTAT(RANK,R,DEV,M,NT,T) 

INPUT PARAMETERS ; RANK, M, NT 
OUTPUT PARAMETERS ; R, T 

REFERENCED BY : The TSTAT subroutine is referenced by the driver routine MAIN. 

SUBPROGRAMS REFFRENCEP ; None. 

INPUT/OUTPUT DEVICES : None. 

LOCAL VARIABLES ; 

Name Type Dime nsion Description 

RANK R 

R R 

Z R 

DEV R 

S R 

T R 


M X NT Array of the ranks of the experimental data. 
1 X NT Rank totals for treatments. 

1 Expected value of R(1,J), J = 1, •••, NT. 

NT [R(1,J) - if- 

1 J CR(1,J) - Z]2 

J 

1 Friedman's T-statistic. 
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M 


ORIGINAL PAGE IS 
OF POOR QUALITY 

3.5 OUTPUT 

The OUTPUT subroutine prints an arbitrary M x N array in the following format: 


"1.1 

XI, 2 

• • • 

XI, 8 

X2,l 

• 

X2,2 

• • t 

X2,8 

* 

Xm,l 

Xm,2 

• • • 

Xm,8 

XI, 9 

XI, 10 

• • • 

XI, 16 

X2,9 

• 

><2,10 

• • • 

X2,16 

• 

"m,9 

• 

Xm,10 

• • • 

Xm,16 

• 

xi,n-q 

• • • 

XI, n 


X2,n-q 

• 

• • • 

X2,n 


• 

Xm,n-q 

• • • 

Xm,n 



where q = M0D(NT,8) - 1. 

CALLING PROCEDURE : CALL 0UTPUT(S,N0,NV) 

INPUT PARAMETERS ; S, NO, NV 
OUTPUT PARAMETERS ; None. 

REFERENCED BY : The OUTPUT subroutine is referenced by driver routine MAIN. 

SUBPROGRAMS REFERENCED : None. 

INPUT/OUTPUT DEVICES: UNIT 6, disk with LRECL = 133, BLKSIZE = 133 
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LOCAL VARIABLES; 


Name Type Dimension Description 


S R 

NO I 

NV I 

H I 

M I 

K I 

L I 


NO X NV Array to be printed. 

1 Number of rows. 

1 Number of columns. 

1 NV 4 8. 

1 M0D(NV,8). 

1 Lower limit of each print iteration 

1 Upper limit of each print iteration 


4. EXAMPLE 


Assuming the FRIEDMAN program has been compiled, the following EXEC routine 
may be used to run FRIEDMAN: 

&C0NTR0L OFF 

FILEDEF 5 DISK &1 &2 &3 (PERM 
FILEDEF 6 DISK FRIEDMAN OUTPUT D 
(PERM LRECL 133 BLKSIZE 13'5 
FILEDEF 16 TERM (PERM 
LOAD (CLEAR NOMAP 
START 

SPOOL E NOH 

PR FRIEDMAN OUTPUT D (CC 

SPOOL E HOLD 

&EXIT 

Suppose the data set to be used is as shown in figure 2, When the above EXEC 
routine is executed, the following prompt will appear: 

ENTER NUMBER OF CARDS PER BLOCK IN 12 FORMAT. 

The user will respond with "SPACE 2." The space is necessary because of the 
12 format. Upon carriage return, the prompt 

ENTER NUMBER OF TREATMENTS IN 12 FORMAT 

will appear. The user should respond with "SPACE 7." Upon carriage return, 
the final prompt will appear: 

ENTER INPUT FORMAT. 

The user should respond with 

(4F3.0,/,3F3.0). 

The output will be spooled to the printer no-hold and should appear as shown 
in figure 3. 
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original page vs 

OF POOR QUAUTY 


14.17.19.14. ) 

Block 1 

10.21. 9. ) 


13.19.15.14. ) 

[ Block 2 
13.13.11. ) 


21. 9.19.17. ) 

} Block 3 

9. 7.19. 


Figure 2.- Data set example. 
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Figure 3.- Output example 



5. REFERENCES 


1. Conover, W. J.: Practical Nonparametric Statistics. John Wiley & Sons, 

Inc., New York (New York), 1971. 

2. Merchant, Michael J.: Applied FORTRAN Programming With Standard FORTRAN, 

WATFOR, WATFIV, and Structured WATFIV. Wadsworth Publishing Company, 
Inc., Belmont, California, 1971, p. 356. 


5-1 


APPENDIX 

FRIEDMAN FORTRAN CODE 


onoon ooono ooonr; onnoo nonno onnon ononn r>r;ono nno 


ORiGPAL FAC^i i;J 
OF POOR QUAUTY 


fill: FMIEOMAM FORTRAN A EODL / JOHNSON SPACE CENTER 


FkIEDmaN'S test for SEVERAL RELAFEO SAMPLE'S 


DIMENSION FMT(80) 
REAL XINPUT(SU0»30) » 
REAL l)EVfl.30> 
iNTEbER INDEX (30) 


KANK(S00»30) . R(l»30) »WURK(30) 


SET UP PROMPTS AT TERMINAL 


WO 

wS 

100 

105 

110 

lib 


format! ‘ENTER NUMHER OF CARDS PER BLOCK IN 12 FORMAT*) 
REA1H16.9S) NC 
FORMAT! 12) 

WRITE! 16»100) 

FORMAT! ‘£NT£R NUMttER OF TREATMENTS IN 12 FORMAT*) 
RFAO!16.10S) NT 
FORMAT! 1 2) 

KVITEI16»110) 

FORMAT! ‘ENTER INPUT FORMAT*) 

READ!16,U5) FmJ 
FORMAT!80Al) 


COUNT NUMBER OF CARD IMAGES DETERMINE NUMBER OF BLOCKS 


M=0 

120 REA0!5tl2S,£ND=5) 
12b P0RMAT!1X) 

M=M+1 
GO TO IBO 
5 M=M/NC 


READ DATA 


CALL INPUT !XlNPUTtFMT»M»NT) 
WRITE!6*130) 

130 format !1H1,‘ input DATA*) 


white input data 


CALL OUTPUT!XINPUT»M,NT) 


SORT TREATMENTS S, ASSIGN RANKS 


DO 20 1=1, M 

CALL SORT!XINPUT,lNUEX,I,M,NT) 

CALL XRflNK !X INPUT tRANK, INDEX, 1 ,M,NT,wORK) 
20 CONTINUE 

WRITE!b,135) 

135 FORMATIlHl,* MATRIX OF RANKS*) 


WRITE MATRIX OF RANKS 


CALL OUTPUT !RANK,M, NT) 
WR1TE!6,140) 

lAO FORMAT!1HO,///» * rank TOTALS*) 
CALCULATE FRIEOMAN‘5 T STATISTIC 
CALL TSTAT !RANK,R»DEV,M,NT,T) 
WRITE RANK TOTALS 


A -1 



noono nnnnn onnon norann oonnn onooo noono 


OF POOR 


IS 

qUAUTY 


FILE: FRIEUMAN FOHTKAN A EuDL 


CALL 0UTPUT(R,1»NT) 

WRITE (eUA'iJ 

145 FORMATdHl,* FWIEUMAN T*) 


white FRIEDMAN'S T STATISTIC 


/ JOHNSON SPACE CENTER 


WRITE(6,150) r 

ISO FORMAT(lHO*//t* T = *fF7.3) 
STOP 
END 


SUBROUTINE INPUT (XlNPUT.FMT»M,NT) 


SUBROUTINE INPUT (XINPUT ,FMT »M«NT ) 
REAL XINPUT (M,NT) »FMT (dO) 

REWIND 5 
l)U 10 1 = 1 »M 

10 PEAD(5»FMT) (X INPUT ( I ♦ J) ♦ J=1 »NT ) 
RETURN 
END 


SUBROUTINE SORT (XINPUT* iNUEXt I tMfNT) 


SUBROUTINE SORT (XINPUT ♦ INDEX* 1 *M*NT ) 
UIMENSION XINPUT (M*NT) ♦ INDEX(NT) 


Initialize index 


JO 10 J=1*NT 
lU INDEX (J)=J 

ndelta=nt 


BtblN SORT 


2U NDELTA=nOELTA/2 

IF (NOElTa ,LT. 1) GO TO SO 
JO NSWlTCsd 

max=nt-ndelta 

DO 40 J=1*MAX 

IF <XINPU7(I*J) .LE. XINPUT(I*J^NDELTA) ) 60 TO 40 
TEMP=XINPUT(I*J) 

ITEMP=1NDEX(J) 

XINPUT ( I . J) =XINPUT ( I » J*N0ELTa) 

INDEX (J) = index (J+NUELTA) 

XINPUT ( I * J*NDELTA) =TEMP 
INDEX (J*NUELTA)=ITEMP 
NSWITC=1 
40 CONTINUE 

IF (N5W1TC .EQ. 0) GO TO 20 
GO TO 30 
50 PETURM 
END 


SUBROUTINE XRanK ( X INPUT *RANK . INDEX, I ,M,NT * WORK) 


SUBROUTINE XPANK (X INPUT *RANK* INDEX* 1 *M* NT* WORK) 
REAL XINPUT (M, NT) , HANK(M*NT)* WORK (NT) 

INTEGER INDEX (NT) 


ASSIGN INITIAL HANKS L STORE IN ARRAY WORK 


00 *0 J=1*NT 
10 WORK(J)=FLOAT(J) 
N=NT-1 
NE0UAL=1 
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OF POOR QUALITY 


FILEJ F^^IEC'MAN FOHTKAN A EODL / JOHNSON SPACE CENTER 


C 

C 

C 

C 

C 


CHECK FOR TIES & FIND MEAN RANK WHEN DETECTED 
DO <20 KsltN 

IF tXINPUTdtK) .NE. XINPUTU»K*in GO TO 30 

NEQUAL=NEQUAL^l 

GO TO 2n 

30 IF (NEQUAL .GT. 1) GO TO AO 
GO TO 20 

AO NFlRSTaK-NEOUAL+1 
INC=NtQIJAL-l 
last*nfirst*inc 

SUM*0. 

, 00.50 E=NFIRST*LAbT 

bO SUM=SIJM*WORK<L) 

DO 60 LU-NFIRST*LAST 
oO WORK (LL ) =SUM/NEOUAL 


C 

C 

C 

c 

c 




20 CONT 


ASSIGN FINAL RANKS TO TREATMENTS 


DO fO MM=ltNT 

70 RANKdtlNOEX (MM) ) =wORK(MM) 
RETURN 
END 


C 

C 

C 

C 

C 


c 

c 

c 

c 

c 


SUHROUTINF OUTPUT (S.NO<NV) 


SUBROUTINE OUTPUT (S. NO. NV> 

INTEGER H 
REAL S(NO.NV) 

H-NV/fl 

IF (H .EQ. 0) Hsl 
MaHOI)(NV.«) 

K = l. 

DO 10 KK=lfM 
, WRITE(6.l05) , 
lob FORHATdX.//) 

L=K*MINO(NV-l.r) 

WRITE(6.110) (N.N=K.L) 
d,0 FORMATdX.3k.a(8X.I2.SX) ) 

WRITE(6.120) 

120 FORMATdX./) 

00 30 1=1, NO 

30 WRITE(6,130) I ♦ (S d . J1 . J=K ,L) 

130 FORMATdHO.I3.H(3X.Fl2,b) ) 

K=L+l 

10 CONTINUE 

IF (M ,EO. 0 ,UR. NV .LE. 8) GO TO 60 
L=K-t«M-l 

WRITE(6,1A0) d.I=K,L) 
lAO FORMATdH0.///»AX.«{HX. I2.5X) ) 
WRITE(6,150) 

IbO FORMATdX,/) 

DO 50 1=1. NO 

50 WRITE{6,130) I . (S (I . J) ♦ J=K.L) 

60 Rf JURN 

SUbROUTINE TSTAT(RANK»R.0£V,M,NT,T>, 


SUBROUTINE TSTAT (RANK.R.DEV.M.NT.T) 
REAL RANK(M.NT) .Rd.NT) 

REAL DEV (1. NT) 


FIND RANK TOTALS FOR EACH TREATMENT 
1)0 30 J=1,NT 
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oonrsn ooooo onnoo 


ORIGINAL Pi\Q.Z IS 

OF POOR QUALITY 


FILh,S pi^IEOMAN FO^^TPAN A tODL / JOHNSON SPACE CENTER 


00 AO I*l»M 

,R(1*J)sW<UJ)*RANK(I»J) 

30 CONTINUt 

CALCULATE EXPECTATION OF RANK TOTALS 
Z*PL0AT{M#(NT*1) >/2. 

CALCULATE SQUahED DEVIATIONS FROM MEAN 


00 so LL*l«NT 

bO DEV(l»LU = (R(lfLL)-Z)»*2 
S=0, 

DO 60 iJsliNT 
60 SsS^DEVn.IJ) 

U=FLOAT(M«NT*(NT*n ' 


CALCULATE T 


T=(I2./U)«S 

RETUKN 

END 


*U.S. OOVERNMENT PfllNTINO OFFICE. lQ«J-5a9-011»« 
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